Fixed points, stability and intermittency in a shell model for 
advection of passive scalars 
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We investigate the fixed points of a shell model for the turbulent advection of passive scalars 
introduced in jjj]. The passive scalar field is driven by the velocity field of the popular GOY shell 
model. The scaling behavior of the static solutions is found to differ significantly from Obukhov- 

1/3 

Corrsin scaling 9 n ~ k n which is only recovered in the limit where the diffusivity vanishes, D — > 0. 
From the eigenvalue spectrum we show that any perturbation in the scalar will always damp out, 
i.e. the eigenvalues of the scalar are negative and are decoupled from the eigenvalues of the velocity. 
Furthermore we estimate Lyapunov exponents and the intermittency parameters using a definition 
proposed by Benzi et al. M. The full model is as chaotic as the GOY model, measured by the 
maximal Lyapunov exponent, but is more intermittent. 

PACS numbers: 



I. INTRODUCTION 

In recent years much attention has been paid to the 
origin of intermittency in fully developed turbulence. 
Progress has been made in the somewhat simpler prob- 
lem of the anomalous scaling in the passive advection 
of a scalar quantity (e.g. temperature or the density of 
a pollutant). Fundamental analytical results have been 
obtained when the advecting velocity field was assumed 
to be Gaussian and delta-correlated in time, the so-called 
Kraichnan model , and in the context of shell mod- 
els Q. Here we consider the perhaps more realistic situ- 
ation where the passive scalar shell model is driven by a 
non-Gaussian velocity field with finite correlation time, 
namely the velocity field of the GOY model In 
this case the model can be analyzed using standard tech- 
niques of (low) dimensional dynamical systems, e.g. the 
study of bifurcations, eigenvalue spectra and Lyapunov 
exponents. 

In absence of convective effects the passive scalar field 
O is governed by the equation: 



d t Q + (v ■ V)6 = L>V 2 9 + F e 



(1) 



Here v is the velocity field, Fq is an external forcing and 
D is the diffusion coefficient. 

According to the analogue of the K41 theory ||] for the 
passive scalar, developed by Obukhov and Corrsin [ p^[ , 
the structure functions 



S p (r) = (\Q(x + r) - 6(x)\ p ) ~ I 



(2) 



(where I — \r\) scale linearly with p, more precisely 
H{p) = p/3. Experimentally, however, one observes for 
p > 3 |ll| strong deviations from this. This is usually 
referred to as anomalous scaling or intermittency. The 
deviations seem to be even more pronounced than for the 



structure functions of the velocity: the passive scalar is 
said to be more intermittent. 

This paper is organized as follows: after having intro- 
duced the model, we examine in section III the scaling of 
its fixed points and in section IV their stability. These 
studies have already been performed for the GOY model 
to which the passive scalar model is coupled. We 
will review these results for the sake of completeness. In 
section V we study the full dynamics of the model and 
investigate its chaotic and intermittent behavior. 



II. SHELL MODEL FOR PASSIVE SCALARS 

Shell models appear to capture many properties of fully 
developed turbulent flows but are easier to study than 
the Navier-Stokes equations (see |lf| for a review). In 
this letter we study the passive scalar shell model pro- 
posed in Ref. The multiscaling of the model is in 
good agreement with experimental data flllfl. The GOY 
model has been studied intensively 
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sive scalar model has attracted somewhat less attention 

Both models are constructed in Fourier space, retain- 
ing only the complex modes u n and 9 n as a representa- 
tive of all modes in the shell of wave number k between 
k n = koA n and k n +\. One uses the following assump- 
tions: (i) the dissipation resp. diffusion is represented by 
a linear term of the form: — vk~Un resp. —Dk^O n (ii) 
the nonlinear terms of the form k n u n >u n " resp. k n 9 n 'U n " 
with (iii) n' and n" among the nearest and next-nearest 
neighbors of n and (iv) in absence of forcing and damp- 
ing conservation of volume in phase space and conser- 
vation of J2 n \ u n\ 2 an d \@n\ 2 - Moreover the scaling 
laws u n ~ k n and 9 n ~ k n ly/3 form a fixed point of 
the inviscid unforced equations. 
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The resulting equations arc llj : 
(— + vk n ) u n = i (a n k n u n+1 u n+2 + b n k n -iu n _ 1 u n+1 

+ C„fc„_ 2 W* l -iU r *_ 2 ) + f$n,4, (3) 

(j t + Dk 2 n ) d n = i (e n «_^; +1 - u* n+1 6* n _ x ) 

+ 5n«-2^-l+<-lC-2) + 

)) + /<W (4) 

A possible choice for the coefficients is: 

a„ = 1 b n = -6 c n = -(1 - 6) (5) 
k n k n —i k n +i 

e n = y g n = — -- hn = ~2~ d 

The boundary conditions are: 

h = b N = ci = c 2 = ajv-i = aw = (7) 
ei = e N = ,9i =g 2 = h N -i = h N = (8) 

For the parameters, we choose for example the follow- 
ing standard values: 

N = 19, A = 2, k = A~ 4 , i/ = 10~ 6 , 
/ = / = 5-10- 3 -(l + z),i?=10- 6 . (9) 

The free parameter S is related to a second quadratic in- 
variant which for the canonical value S = \ is similar to 
the helicity |§. 

These equations determine the evolution of the vector 
(U, O) = (Kui, . . . , SRwat, 3uat, K^i, . . . , 36*at) and 
thus form a 4iV dynamical system. 

III. SCALING OF FIXED POINTS 

A first step towards a full understanding of the model 
consists of an investigation of its static properties. In 
this section we examine the scaling properties of the fixed 
points of equations (j^)-fl) ■ The major problem is to find 
the static solutions of (||), those of (Q) can then be found 
easily because 8 is linear in both u and 9. 

It was found in |l^] that (||) in the unforced inviscid 
limit allows three self-similar static solutions: the trivial 
fixed point u n = 0, a "Kolmogorov" fixed point u n = 

— 1/3 

k n gi(n) and a "flux-less" fixed point u n = k~ z g2(n), 
with gi(n) and g2(n) any function of period three in 
n and z = (— lii\(S — 1) + l)/3. The corresponding 

fixed points of (Q) are: 9 n = 0, 9 n = kn 1 ^ 3 gi{n) and 

9 n = k n 2 . We will here focus on the Kolmogorov 
fixed point which is believed to be the most important 
for the dynamics ]l9| although it was suggested that also 
the trivial fixed point might play a major role pT| . 

We note that the static solution for u n = u n e^" can 
be turned into real form by a change of phase. Following 
Schorghofer et al. [|l4| we choose the phases: 

C \tx for n = 1,4,7,... 
4>n = < Itt for n = 2,5,8,... (10) 
[ |7r for n = 3, 6, 9, . . . 



It can be shown that the static solution of 9 picks up the 
same phase as that of u. 

As has been observed ]l9|] , the dynamics of the sys- 
tem converges to the Kolmogorov fixed point for S < 
0.379634. When increasing 6 the system undergoes a 
series of Hopf bifurcations and becomes chaotic at S = 
0.38704 (l|Jl|Q. In order to find the Kolmogorov fixed 
point one can thus vary d in small steps and refine the 
solution with Newton's method 

To study the scaling behavior of the static solutions, 
we apply a much larger number of shells, using the same 
parameter values as Kadanoff et al. f fuj : 

N = 90, A = 2, fc = A -1 , v = D = 10" 3 16" 26 , / = / = 1 

(11) 

The forcing acts in this case on the first shell. 

For the static solution of u n we obtain the same result 
as Q, see figure ||, where \og 2 {u n+ \ / u n ) is plotted ver- 
sus n. It is averaged over a period three to get rid of the 
well known period three oscillations. The solution is seen 
to follow Kolmogorov scaling. 




FIG. 1. Static solution of the GOY model at S = 0.5 
for the parameters (11). To improve the scaling behavior, 
log 2 (u„+i/u n ) is averaged over three consecutive shells. 

Surprisingly however, the fixed point for 9 n (with the 
fixed point of u n inserted) deviates strongly from the 
Obukhov scaling for a finite value of the diffusivity D, 
as can be seen in figure |[ There is clearly not a well- 
defined power law scaling. This solution also displays the 
period three behavior. In the diffusive range the solution 
looks somewhat noisy, presumably due to the boundary 
conditions. It is clear, that there is a "slow" bending in 
the diffusive regime but as D approaches zero, the curve 
becomes more and more flat and we recover the Obukhov 
scaling in the limit D —> 0. We thus note that in contrast 
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to the velocity-case where the viscosity only affects the 
viscous range, for the passive scalar the diffusion seems 
to act on the whole inertial range, at least for Prandtl 
numbers Pr = -jy ~ 1. This might have its origin in the 
linear character of the problem. 



The matrices D nm and C nm are the contributions of the 
dissipation and cascade terms, respectively. Their ex- 
pressions are: 



A. 



i / k n S nm 



(15) 



and 



-0.2 




FIG. 2. Static solution, fixed point, of the passive 
scalar model where log 2 {0 n +i /&n) (averaged three consecutive 
shells) is plotted versus the shell index n. Here, the number 
of shells is iV=61. 



IV. SCALING OF EIGENVALUE SPECTRA 



du, 



-(k n u n+1 u n+2 - 8k 



+ 1 



-(1-S)k 



(16) 



where the index (0) has been dropped, for convenience. 
The linearized equation for the variation of the phase 

is: 



(17) 



Thus the stability eigenvalues of modulus and phases 
differ only in a minus-sign in front of the cascade term. 
We obtain similar results as in |14| ). For values of 8 < 
Sbif « 0.37 all the eigenvalues of both phase and modulus 
matrix have negative real part. Above this value, some 
of the real eigenvalues of the phase matrix turn complex 
(in pairs) and cross the imaginary axis. At 8 = 0.5 they 
have become real again, but now positive. This situa- 
tion is shown in figure || The eigenvalues of the modulus 
matrix eventually turn positive at 8 « 0.7. In figure ^ 
we have multiplied both real and imaginary part with a 
factor p 0: 



P 



log 2 (l + 2 10 M) 



(18) 



Thus the phase remains unchanged, while the modulus 
is rescaled. In case the eigenvalues have constant ratios, 
they are evenly spaced on the plot. Thus we are able to 
visualize both the very small and very large eigenvalues. 



Now the stability of the Kolmogorov-Obukhov fixed 
point is examined in terms of the eigenvalue spectra. The 
system eq.&El) is written as 



(12) 



u = f(u) 
6 = g{u,6) 

The Jacobian matrix is symbolically given by: 



J 



du ( 
dg{u,e) 

du 




(13) 



The matrix ^l"' 9 - 1 will not matter for the eigenvalues 

cm ° 

of J. The eigenvalues of were studied in (l4|] and 

p3| . It is most convenient to look at disturbances of 
phase and modulus of the velocity variable u n = 



a i<t>n 



,(°) 



l(0) 



u n = Un' + 8u n and </>„ 
for the stability of the modulus: 



One then obtains 



8u r: 



(14) 
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FIG. 3. Eigenvalues of phase matrix (O) and amplitude 
matrix (□) for the velocity field in a polar plot for 8 = 0.5. 
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It is quite trivial to generalize this method to the pas- 

d9 



sive scalar in order to calculate the eigenvalues of dg ^ u ' e ^ 



and one finds: 



X D nm + Bnm) 89 „ 



(19) 



where Df lm is the dissipation term, very similar to the 
one before 



Tr — Dh^-ft 

^nm ~ dJn 'n u 'nm 



(20) 



and B nm is the contribution of the cascade term given 
by: 

1 d 

Bnm = ~ „„ (k n {u n -l9n+l ~ Wn+l^n-l) 
Z OO m 

— k n -i(u n -20 n ~l + U n -l9 n -2) 
+ 2 fc "+l( M n+2#n+l + U n+ i9 n+2 )) (21) 



We see that B nm does not depend on 9 as the 9- variation 
will be differentiated out. 



For the phase disturbance 8ip n one gets 



(22) 



We note that since B„,„+2 = ^K+lUn+i = -B n+2 , n 
and B n ,n+x = \{k n +\u n +2 + = —Bn+i, n , the 

matrix B nm is antisymmetric and its eigenvalues will 
be purely imaginary. The stability of phase and mod- 
ulus will thus be the same. This is actually a conse- 
quence of the conservation of 1 6*^ | 2 . In the model any 
cascade term a n u n + n >9 n +n" has to be supplemented by 



a term 



'U n + n i - n "9 n -n" , since the conservation 
= 0. The first term gives a contribu- 



implies ^n#? 
tion to the matrix B nm : B n , n +n" = a„u„ + „<, the sec- 
ond B n<n - n ri = — a n - n 'iU n + n > - n » which corresponds to 
B n +n",n — ~ci n u n+n i . The matrix is thus antisymmetric 
and this is not an artefact of the model since it stems 
from a conservation law also valid in a real system. 

One does not expect to find any bifurcations since the 
diffusion matrix —D nm will cause the eigenvalues to have 
negative real parts, whatever the value of the driving ve- 
locity. 

Indeed one finds that for all values of 8, the spectrum 
of eigenvalues of the phase matrix is like shown in figure 
[|, where the eigenvalues are again multiplied by p. One 
observes evenly spaced eigenvalues organized in branches. 
The presence of three branches (in both upper and lower 
half of the complex plane) might be caused by the period 
three of u n (in n). If one inserts in the matrix B nm a 
fixed point of an imaginary model without period 3, one 
just finds one branch in each half plane. 
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FIG. 4. Eigenvalues of phase matrix (O) and amplitude 
matrix (□) for the passive field in a polar plot at S = 0.5. 

Let's us now consider what the above imply for the 
dynamics of the model. At time t — 0, starting from 
the particular state of the scalar 0„(O) we impose a small 
perturbation: 9' n (0) — 9 n (0) + 89 n (0) and look how the 
perturbation evolves in time, subjected to the same ve- 
locity field. For modulus and phase of 89 we have again 
equations ( [l9| ) and (22), since the matrices D 9 and B 
do not depend on 9. Because the eigenvalues of B — D 
have negative real parts, the perturbation will damp out, 
meaning that after some time 9' n (t) ~ 9 n (t) . This has 
been observed already by Crisanti et al. 



V. INTERMITTENCY 



The term intermittency is used in different contexts 
and a precise mathematical definition does not exist. In 
turbulence one speaks of intermittency corrections to the 
Kolmogorov or Obukhov-Corrsin power law. In dynami- 
cal systems in general, intermittency means the presence 
of quiescent periods randomly interrupted by burst. It is 
believed that these phenomena are related: the scaling 
corrections should have their origin in intermittent be- 
havior (in space and/or time) of velocity or energy dissi- 
pation. Here we ask ourselves the question whether the 
more pronounced deviations from classical scaling for the 
passive scalar are reflected by a more intermittent behav- 
ior (in time) of the passive scalar field. In order to test 
this we invoke a definition of intermittency for dynamical 
systems proposed by Benzi et al. 

Firstly, the response function i?t(r) is defined as the 
rate at which a disturbance vector Sx(t) of the system 
x = f(x) has expanded after a time r. 



Rt{r) 



\8x{t + r)\ 
\8x{t)\ 



(23) 
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In chaotic systems one typically observes that 
\6x(t + t)\ ~ \5x(t)\e rX , where A is the maximum Lya- 
punov exponent. This maximum Lyaponov exponent can 
be calculated by averaging the instantaneous Lyapunov 
exponents lni M r L \ — ]irti r _ too QsJMlll , 

Intermittency can be thought of as connected to the 
fluctuations of the instantaneous Lyapunov exponent. 
An intermittency parameter fi is thus defined as the vari- 
ance of In R(t) |2|: 



<(lni?(T)) 2 )- (In R(t)) 2 =llt 



(24) 



Numerically we proceed as follows: the equations of 
motion are numerically integrated along the equa- 

tions for the disturbance vector Sx = J{x)8x where 
(x) = (u, 9) and J is the Jacobian matrix. We then esti- 
mate the response function and the first two cumulants 
of its logarithm. From time to time the disturbance vec- 
tor is renormalized in order to avoid numerical overflow 

Consider now the Lyapunov exponents of the full sys- 
tem of velocity and passive scalar field ([|^|), with the 
parameters (^|). We distinguish how the perturbation, 
initially applied equally on both systems (8u n = S9 n = 
constant for all n), evolves independently in each of the 
two parts of the system. Therefore we introduce response 
functions that measure the expansions on the velocity 
part r[ u \t) and on the scalar part rI S \t). They are 
defined as: 



(„) _ \Su(t + T)\ 

4 {T) -^(tr 

{e) \S9(t + r)\ 

1 {T) -^mr 



(25) 



The corresponding maximal Lyapunov exponents is de- 
fined as: 



x(u) = (lnflW(r)) 



(26) 



and the analogue for \( e > . This is reminiscent to the "Eu- 
lerian Lyapunov exponent" and "Lagrangian Lyapunov 
exponent" introduced by Crisanti et al. | p6| , p3[ . However 
the Lagrangian behavior of the particles is only equiv- 
alent to the Eulerian passive scalar field in the case of 
vanishing diffusion. 

The intermittency parameters are related to the vari- 
ance of In R(t) is the same way as before. We find nu- 
merically: 



X (u+e) _ A («) w A (9) = 0J65 ± _ 002 



(27) 



where we have denoted the Lyapunov exponent for the 
full system as: \( u+e \ On expects theoretically: 



A («+e) 



t[A< u >,AW] 



(28) 



since the disturbance vector will evolve towards the most 
expanding direction. Our result is obviously in agreement 
with this relation. For the Eulerian and Lagrangian Lya- 
punov exponent, one has numerically found the generic 



inequality Xl > Xe- Our result X^ w X^ might be due 
to the presence of a finite value of diffusivity. 

The intermittency parameters /i differ on the other 
hand significantly from each other. We find (for S = 
0.5): = 0.127 ± 0.003, n {u) = 0.125 ± 0.003 and 

p^ 6 < — 0.151 ± 0.003. Thus the passive scalar behaves 
equally chaotic, but more intermittent, than the veloc- 
ity. In order to understand this we write the equation for 
9 as follows: 



= A(u)6 + f 



(29) 



where A is a matrix depending on it, whose eigenvalues 
a fulfil the inequality Jier < 0. Thus, if u is constant in 
time, 9 will converge to its fixed point 9 = — A(it)" 1 /. 
Alternatively, if u fluctuates as for S = 0.5 in the GOY 
model, this "fixed point" will fluctuate as well. More- 
over the convergence towards it will be irregular since 
the eigenvalues of A(u) vary. It is therefore natural to 
expect that the behavior of 9 is more intermittent than 
that of it. 



VI. CONCLUSIONS AND OUTLOOK 

We studied various properties of a shell model for the 
advection of a passive scalar. We calculate the eigenvalue 
spectrum of the Kolmogorov fixed point and show that 
the passive scalar part is stable against perturbations. 
This is in fact not very surprising since it follows from 
conservation of @ 1 2 ■ 

The fixed point does on the other hand not follow the 
Obukhov-Corrsin scaling. This suggests that the scaling 
properties of the passive scalar are much more sensitive, 
than the velocity, to non-inertial properties such as the 
dissipation and perhaps also the forcing. In experiments 
the scaling zones of the passive scalar are less clearly ob- 
served and some "nonuniversality" (with respect to forc- 
ing) seems to have been observed in [ p7| , but we of course 
cannot pretend that the mechanism observed here will be 
found also in real systems. 

We have measured the relation between Lyapunov 
exponents and intermittency parameters of the passive 
scalar and velocity. We find: A^"^ rj X^ and 



(0) 



Thus the passive scalar behaves more intermit- 



tently which is in agreement with its more pronounced 
deviations from linear scaling. It would be interesting to 
establish these relations from the equations of motion. 

It could also be interesting to investigate the relation 
between coherent structures in the GOY model [^1 20 
and those in the passive scalar model. 
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